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■ Abstract 

We consider Majorana neutrino conversions inside supernovae by taking into account both 
' flavor mixing and the neutrino magnetic moment. We study the adiabaticity of various possible 

■ transitions between the neutrino states for both normal and inverted hierarchy within the various 
solar neutrino problem solutions. From the final mass spectrum within diffrent scenarios, we 

I , infer the consequences of the various conversion effects on the neutronization peak, the nature of 

final spectra, and the possible Earth matter effect on the final fluxes. This enable us to check 
£^ ' the sensibility of the SN neutrino flux on magnetic moment interaction, and narrow down possible 

scenarios which depend on: the mass spectrum normal or inverted, the solution of the solar neutrino 
• problem; and the value of /iB. 

m : 
o ■ 

1 Introduction 

The neutrino signal detected on the Earth from the SN1987A explosion T; has opened up new ways 

to probe the neutrino properties. Initially only constraints on the static properties of neutrino during 

its interstellar journey (namely from vacuum oscillations) could be acquired. Now with the entry into 

service of the large neutrino underground detectors like Superkamiokande (SK) j2j and the Sudbury 

X ' Neutrino Observatory (SNO) 0, modification of the flavor content of the neutrinos emerging from 

\—{ ' 

the SN itself could be contemplated. This is helped greatly by the new quantitative progress achieved 
in the past few years in our knowledge of the mass squared differences and mixing parameters from 
the atmospheric neutrino oscillations (SK and MACRO [I] in particular), and the solar neutrino 
oscillations 

Various works have appeared studying the consequences of both the matter resonance effect (the so 
called MSW effect) (Hj, and the spin flavor precession (the RSFP effect) (Zj, especially with the present 
data on the neutrino magnetic moment jSj, due to the very large magnetic field in the pre-supernova 
interior These studies have been mostly confined to the case of two active neutrinos. The case of 
sterile neutrinos is no more rigorously pursued in the light of the strong constraints on their existence 
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from the atmospheric neutrinos experiments. They used to be invoked in neutrino conversion, in 
particular for their potential role in enabling r-process nucleosynthesis jlUj . 

We wish in this paper to consider the general case of three active neutrinos on the neutrino SN 
spectrum, namely the ones associated with the three known leptonic flavors, in the context of both 
the MSW effect and the RSFP effect. We have followed the seminal work of Dighe and Smirnov 
as far as their general line of attack on the question, but generalizing it to take into account both 
of the above mentioned effects at the same time. Considering the various schemes of neutrino mass 
and mixing allowed by the data, we systematize the discussion of the various spectrum distortion 
effects, and thus check the discriminating power of such studies and how they could help resolve the 
ambiguities associated with it, notably: 

• the solar neutrino solution (which one to choose?) 1 . 

• the type of hierarchy for the neutrino masses {normal or inverted). 

• the value of fiB± where /i is the neutrino magnetic moment, and B± the SN transverse magnetic 



It follows from the existing solar and atmospheric neutrino data that neutrino mass-squared differences 
satisfy the hierarchy Am^ <C Am^, which permits us to order the three mass-squared according to 



be negative (See Fig.l-b). We will assume that Am^ is relevant for the solar neutrinos oscillations, 
while Ato^ is relevant for the oscillations of the atmospheric neutrinos. 

The effects of the neutrino conversions can be observed through, (i) the disappearance (partial or 
complete) of the neutronization peak; (ii) the interchange of original spectra and the appearance of 
a hard v e spectrum; (Hi) the modification of the v e spectrum; (iv) the Earth matter effect, which is 
studied taking into account neutrino mass and mixing in I lj. Note that for significantly cosmological 
mass-squared differences (Am 2 = I ~ 100 eV 2 ), the spin- flavor precession and resonant spin- flavor 
conversions may affect the supernova shock reheating and r-process nucleosynthesis (U]. 
Although the ambiguities on the neutrino spectrum could not be solved, a systematic study of their 
effects in the general three-generation case pave the way for further constraining the various scenarios 
as new limits are obtained from the detectors on the Earth. 

This paper is organized as follows: we first obtain the original neutrino fluxes emerging from the 

SN core; we then study neutrino conversion outside the core for both normal and inverted mass 

hierarchy, and finally determinate the final neutrino fluxes which reach the Earth detectors. The 

various neutrino conversion effects are then classified according to the different parameters relevant 

for the Solar neutrino solutions and the magnitude of the magnetic interaction term fiB±. 

x The recent results from the KamLAND experiment |12| indicate that the LOW scenario is most probably ruled out, 
which leaves the LMA as the most favored one. 



field. 
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Fig 1: Neutrino mass pattern: a) represents the case within normal mass hierarchy; and b) represents 
the case within inverted mass hierarchy. 



2 The Original Spectrum 

When the inner iron core of a massive star going through a type // supernova explosion becomes unable 
to support the electron degeneracy pressure, most of the gravitational binding energy is released in 
the form of a violent burst of neutrinos of all species |131 114j . 

These neutrinos produced both during the neutronization burst and the subsequent thermal cooling 
could undergo transformation of kind both inside the SN, and outside on their way to the Earth. We are 
however interested only in their transformation within the SN. In some scenarios these transformations 
could boost the energy deposition at the stalled explosion front as well as play a key role in the explosive 
nucleosynthesis [T^l- Yet, since in our study we are using mass differences relevant to the Solar Neutrino 
Problem (SNP) and the Atmospheric Neutrino Problem only, these neutrinos transformations are not 
expected to change the dynamics of the explosion as the corresponding resonances take place outside 
the core. 

Two effects have been widely considered: the first one is the matter effect, i.e. the interaction of z/s 
with different matter constituents which are mainly electrons, protons and neutrons. The second one 
consists of the interaction of the neutrino's magnetic moment with the transverse magnetic field of 
the SN. For Dirac neutrinos, they flip into a right handed vr kind which is known to be sterile since 
undetected till now. Comparing the inferred energy output of SN1987A to the theoretical expectations, 
one can put strong constraints on these Dirac neutrino conversions. On the other hand, when including 
the most general mass terms in the simplest extension of the Weinberg-Salam model, the neutrino 
fields, once diagonalized, turn out to be of Majorana kind. All this concurs to make the Majorana 
neutrinos more palatable than the Dirac kind in the context of SN studies. We will assume in what 
follows that the neutrinos are Majorana particles. Now since they are their own antiparticles, only 
flavor conversion will be allowed. 

In order to study these effects together, we should identify the profile of both the density and the SN 
magnetic field. We take the standard parametrization: 




(1) 
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where B± is the magnetic field strength at the distance r o =10km; and k=2 or 3; in our work, we will 
consider k = 3. The value of B± lies between 2 (10 12 ~ 10 15 ) G , where G is the strength magnetic field 
unit. For this density profile, the important effects occur between 3 0(lgcm -3 ) and O(10 7 gcm~ 3 ), 
i.e. this region may not be affected by the SN shock wave. We will thus take the progenitor profile 
for p > 1 g cm -3 as approximately given by ^Z] 

P Y ^Po{k)" Z (2) 

where p Q =2x 10 13 g cm -3 . Of course, the exact shape of the density profile will depend on the precise 
composition of the star. 

The first question to answer is the flavor content of this burst when reaching the conversion regions. 
The key point here is that they are produced in a high density medium so that their mass eigenstates 
can readily be inferred. 

For this matter, let us rewrite the Hamiltonian in the flavored basis v* =(iy e ,u^,u T ,i> e ,i> f j i ,i> T ) which is 
related to the mass eigenstates basis in vacuum i^=(z^i, ^2,^3, ^1,^2)^3) by the unitary transformation 
= Uu v . We can then deduce the initial mass spectrum from the evolution of neutrinos at very high 
density (p > 10 s gcm~ 3 ). 

Taking into account the neutrino effective potential coming from both the MSW effect and the RSFP 
effect, we can cast the Hamiltonian into the following form: 

H = Vacuum term + Matter term + Spin Precession term (3) 

where each term, after the substraction of terms proportional to the unity, is given by: 

Vacuum term = 

*" \ U ) \ M z J \ W J 
_ v ^;„„ Jv n n 1 ov 1 V 1 V \ 

(4) 




Matter term = ap x diag {Y e ,0,0,l-2Y e ,l-Y e ,l-Y e } 

( 



Spin Precession term = ( ) , where C 
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with M 2 =diag{0, Am^, Am 2 ^}, a = \f2Gpm^ , Y e = - n ^ n is the electronic fraction, the Am 2 's 
are the difference between squared masses; p e/I is the transition magnetic moment between the two 
flavored states v e and and U is the mixing matrix given by ^B] 

C12C13 S12C13 S13 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 1 ( 5 ) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 



2 This is taken to be the transversal value of the magnetic field strength at the surface of the neutron star. 

3 The resonance densities are approximately given by ~ mN 2^G C °e ^ q' wnere q is Y e for the MSW transitions, and 
1-£Y e for the RSFP ones. The values of Am 2 and 8 are taken from the solar and atmospheric neutrino data. Y e is 
generally set to be half, and therefore 1-2Y e lies between 10' J> to 10 then all transitions occur in a range extending 

from a few g . cm -3 to 10 1 %. cm -3 . 
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where Sij = sin6ij,Cij = cosdij for i,j = 1,2,3 (i < j); and 9ij is the vacuum mixing angle between v$ 
and fj. The mixing matrix terms satisfy the unitarity relations: |L^/| 2 = ^ \ Uu\ 2 = 1. The equality 
in the second term in Eq. ©, however, is only true to leading order since radiative corrections induce 
tiny differences between the neutral current potentials of v e , Vn and v T , and in particular, results in 
a very small v^-v T potential difference | — | ~(1CP 5 ~ 10~ 4 )xap |19| . 

Performing a rotation in the non-electronic subspace as {v^,u T ) — > {v^,Vt') which diagonalizes the 
vacuum term and leaves the matter term invariant, while for the magnetic term, the two values p efl 
and p eT become nearly maximally mixed, and so whatever their values, the rotated ones will be of 
the same magnitude, and therefore we have p efl > ~ p eT > ~ p. As for the third one p^ T it remains 
invariant. Since and u T are both produced via neutral currents only and are indistinguishable as 
far as their detection, we take the non-electronic original fluxes to be equivalent, so we can then write 

Comparing at very high density, the values of the various Hamiltonian terms ; we notice that for 
energies in the range of MeV, ap — 2 T<P>~2"E 1 - Since we took the term p to be between [2*U] 10~ 12 /ig 
and 10~ 10 p B , the same thing applies when we compare the matter term with the magnetic one with 
respect to the chosen range of pB± D i.e. ap 3> pB± Q . Note that the chosen values of both p and B± a 
don't affect the SNP solution. We then deduce that at the SN core, the matter term is the dominant 
one, and therefore the Hamiltonian is approximately diagonal 

H ~apx diag { Y e ,0,0,l-2Y e ,l-Y e ,l-Y e } (6) 

This means that the matter eigenstates coincide with the flavored states, and so, the neutrinos emerg- 
ing from the SN core are as follow: 



1- Normal Mass Hierarchy 




~o _ 
v l ~ V e 














2- Inverted Mass Hierarchy 




~o _ 

v x ~ iv 















3 The Neutrino Flavor Dynamics 

In order to find the final spectrum, we should study all possible transitions between the neutrino 
species. We note that these transitions occur in the isotopically neutral region^, which consists mainly 
of layers of 4 He, 12 C, 16 0, 28 Si and 32 5 whose nucleus have N ~ Z, so that the electronic fraction Y e 
will be close to half. We also consider the following parameters values which are: first, when the SNP 
solution is the LMA scenario, the present parameters are sin 2 29 & = sin 2 2#i2 = 0.91 and Am^ = 
6.9 x 10~ 5 eV 2 , while for the LOW scenario, we take sin 2 26 l2 = 0.92 and Am^ = 1.3 x 10~ 7 eV 2 . 
Secondly, for the atmospheric data, we have sin 26 atm = sin 2 26*23 = 1 and Am| 2 = 2.7 x 10~ 3 eV 2 . 
4 The region characterized by l-2Ye=10~ Jt ^10~ !1 , is called the isotopically neutral region. 
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Thirdly, for the 1-3 mixing, we choose a value just below the experimental limit, which is given by 
CHOOZ and satisfies the adiabaticity condition at the H resonance, which leads us to take 
sin 2 (9i3 = 1(T 3 . 

Neutrino Resonances 

The transition between possible neutrino states occur at some preferred regions, called resonance 
regions, which are characterized by the equality between two Hamiltonian's diagonal elements. Thus 
the possible transitions are of two kinds: the MSW kind occurring due to changing electronic density 
when both diagonal elements corresponds to two neutrino states or two anti-neutrino ones, and the 
RSFP kind corresponding to resonant spin flavor in presence of a magnetic field [7j and involving a 
neutrino and a anti- neutrino state of different flavors. In general we expect 4 significant resonances, 
two of each kind namely: the known two MSW ones which are called L and H transitions, and the 
corresponding RSFP ones which we call L and H transitions; in addition to a fifth one of RSFP kind, 
which we call A, but which doesn't affect the neutrino dynamics. Let us see how we can 
deduce these four resonances from general considerations on the full matrix. Starting from the 6x6 
matrix with its six diagonal elements to be equated two by two, we have Cg=15 combinations. Out of 
these 15 combinations, 3 can be eliminated due to the Majorana character which doesn't allow for the 
like-flavor spin flip, 2 can be further eliminated due to the smallness of the v^-Vr potential difference 
and likewise for u^-i> T . For the remaining 10 cases, they constitute five pairs of conjugate transitions, 
for which there could only occur a transition from each pair at a time; since the resonance density 
computed by equating diagonal elements can be of both signs, while it's conjugate transition will be 
necessarily have the inverse sign, and since physically only the case with positive value is relevant. 
This leaves us with five transitions, the fifth one being as mentioned previously not significant. We 
will study each of the MSW and RSFP kinds separately. 
MSW Resonances 

In general, MSW transition between two neutrinos (or anti-neutrinos) is studied around its resonance 
and the general form of the 2x2 submatrix is: 




(7) 

where z = 2apE, and the resonance corresponds to z res = The transition will be adiabatic, i.e. 
there is no jumping between eigenstates, for 7 res > 1, where 7 res is the adiabaticity parameter [22] at 
the resonance given in this case by: 

lres={ E wY (8) 



with 

E na = 2.74 x 10 9 x y} x (JUL V x (JsUj^MeV (9) 
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The parameters b and c are given by: 



Mass hiera- SNP sol- b (L) c (L) b (H) c (H) 

LMA 3.29 x 10~ 5 1.81 x 10^ 8.46 x 10^ 2.67 x lO' 3 

LOW 6.31 x 10- 8 -2.60 x 10~ 6 8.53 x 10~ 5 2.<JS x 10~ 3 

LMA -3.29 x 10' 5 2.35 x 10- 5 -8.61 x W' 5 -2.72 x lO' 3 

LOW -6.13 x 10~ 8 2.80 x 10~ 6 8.53 x 10~ 5 -2.69 x 10~ 3 



Normal 
Inverted 



b and c are given in eV 2 . The sign (-) for the c value means that the conjugate transition occurs 
instead of the usual one, for example in the case of inverted mass hierarchy within LMA solution, 
c(H)=-2.7x 10~ 3 eV 2 means that v\ < — > i>z occurs instead of v\ < — > v% 5 . 
The values of E na are given in MeV at each layer within all scenarios by: 

Mass hierarchy SNP solution L H 

LMA 2.11 xlO 5 I64.6 



Normal 
Inverted 



LOW 7.20 xlO' 2 166.2 

LMA 1.25 xlO 5 167.2 

LOW 5.69xl0- 2 166.2 



RSFP Resonances 

RSFP transitions occurring between a neutrino state and an anti-neutrino one is also studied around 
its resonance and is of the form: 

lie sz \ 
2.E \ sz dz I 



(10) 



but with s = ^ ~1.9xl0 9 x ( ) • Likewise, we find that the adiabaticity parameter is given as 
in Eq. ©, but with 

E na = 1.89 x 1(T 17 x eH x (J^) x (ng^MeV (11) 

where d =l-2Y e for both L and H, d=l-Y e for A. Note that E na is proportional to (/J-B± ) 3 , which 
shows that the adiabaticity of the RSFP resonances depends strongly on the interaction term fiB± . 
The values of c is given in eV 2 . 

Mass hierarchy SNP solution c (L) c (H) c (A) 

LMA 1.81 x 10~ 5 2.67 x 10~ 3 2.65 x 10~ 3 

LOW -2.60 x 10~ 6 2.69 x 10- 3 2.71 x 10~ 3 

LMA 2.35 x 10~ 5 -2.72 x 10' 3 -2.74 x W' 3 

LOW 2.80 x 10~ 6 -2.69 x 10~ 3 -2.71 x 10~ 3 



5 Note that the L resonance transition occurs in the antineutrino channel for the case of normal mass hierarchy within 
the LOW scheme (For more details see Appendix A). 
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Fig 2: The dependence of the flip probability on the neutrino energy. 



We give in the following table, the values of E na in MeV multiplied by ( ^gf ) 3 for L and H, and by 



( ^ b q° )^ for A, at each layer within all scenarios: 

Mass hierarchy SNP solution L H A 

LMA 3.38X10- 11 4.99xl0- 9 2.83xl0~ 19 

N OTTflCbl 

LOW 4.86xl0~ 12 5.02xl0- g 2.90xl0- 19 

LMA 4.39xl0- u 5.08xl0~ 9 2.93xl()- 19 

LOW 5.23xl0~ 12 5.02xl0~ 9 2.90xl()- 19 

Note that stating that the resonance corresponds to the equality of diagonal elements of Eq. Q is not 
correct strictly speaking. The resonances, which corresponds to the minima of the difference between 
the Hamiltonian eigenvalues, match up exactly the diagonal elements when the off-diagonal elements 
are r— independent, which is the case for the MSW effect. On the other hand, the resonances will be 
shifted when the off-diagonal elements are r-dependent as in the case of the RSFP effect. In our case 
however where the r-dependency is of the form p oc B ~ r~ 3 , this shift is negligeable as we show in 
Appendix B. 

The Adiabaticity at Neutrino Resonances 

At the resonance, jumping probability is given by [22] 

Pf = exp 



{-f[^] f } (12) 



which depends on the energy E as we see from Fig. 2. 



Now the observable part of the supernova neutrino spectrum lies mainly between the energies 5 and 
50 MeV and we will consider only energies in this range in our work. One can then divide the whole 
range of energy in three parts: 

• Part I: Pf < 0.1 ~ 0, corresponding to E na > 70 MeV, and where pure adiabatic conversions occurs. 
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• Part II: 70MeV > E na > 0.2 MeV, for which 0.1 < P f < 0.9. In this range, P f increases strongly 
with the neutrino energy. The adiabaticity is partially broken. 

• Part III: Pf > 0.9 ~ 1, corresponding to E na < 0.2 MeV, where the flip probability is close to 1, and 
which leads to a strong violation of adiabaticity. 

Then one can find easily that the transition at the H resonance is completely adiabatic for all possible 
scenarios, as mentioned above. While, at the L resonance, the transition is: completely adiabatic 
within LMA scenario, and therefore Pl ~ 0; and completely non-adiabatic within LOW one and 
therefore Pl ~ 1, for both normal and inverted hierarchies. 

For RSFP transitions, the adiabaticity depends on the value of fiB± , and therefore we will divide the 
whole range of fJ,B± according to the adiabaticity at the H resonance, then the fJ>B± regions are: 

I- Completely non-adiabatic if: nB± < 340fi B G. 

II- Partially broken if: 340^ B G < fiB ±0 < 2400^ B G. 

III- Purely adiabatic if: nB^ Q > 2400fi B G. 

We have plotted in Fig. 3 the flip probability at the H layer as a function of the energy for different 
values of /j,B ±0 . The curves from (1) to (6) correspond to ^g^=300, 500, 1000, 1500, 2000, 2500 
respect vely. 



The dependance of the flip probability at the layer H, on the neutrino energy for different values of /J,B± Q . 

At the A resonance layer, the E na values are very small within the /j.B± q considered range, which 
corresponds to a complete non-adiabatic conversion, from which we deduce that the transition at the 
layer A doesn't affect the neutrino flavor dynamics. For the L layer, the adiabaticity depends on the 
SNP solution and the specific mass hierarchy. Thus, if fJ.B± Q is larger than a given value (p>B± )2, the 
transition will be purely adiabatic; if it is less than the other limiting value (nB± )i, the transition will 
be completely non-adiabatic; and finally when fiBj_ lies between these two values, the adiabaticity 
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SNP solution 


(/U-B_Lo)l 


0-B_l o ) 2 


LMA 


1800 


12700 


LOW 


3500 


24300 


LMA 


1100 


11600 


LOW 


3400 


23700 



will be partially broken. We give in units of [i B G, these two limiting values for the various scenarios: 

Mass Hierarch 
Normal 

Inverted 

Note that the ^B± interval where the L (77) adiabaticity is partially broken intersects with the 
similar one for the 77 (L) resonance within the LMA solution, but it doesn't within the LOW scheme. 
Then the region 77 can be divided in two parts 6 : the first one corresponds to the completely non- 
adiabatic resonance at the L layer, i.e. fJ-B± lies between 340fi B G and 1800 [i B G (1700 n B G for the 
case inverted hierarchy); while for the remaining part from 1800 fi B G (1700 fi B G) to 2400fi B G, it 
corresponds to the partially broken adiabaticity at the layer L. 

For this purpose, we will give the level crossing for each one of different scenario, which are given 
in Fig's 4, 5, 6 and 7, We will check later where the transition occur, namely whether is it in the 
neutrino or anti-neutrino channel (for MSW), and which one of the two RSFP conjugated transitions 
is involved. 

The hierarchy of the densities of the resonance layers (L with 77, and therefore L with 77 too), leads 
to the factorization of the neutrino flavor dynamics, the transitions between two of resonance layers 
can be considered independently, and therefore each transition is reduced to a two neutrino problem 
jllj . Furthermore, the MSW-RSFP resonances regions non-overlapping condition, for L-L or 77-77 
transitions, is given by |2l] 



L p |tan20 Q/ g| + 



2Va8 B ± ( r i) 



Am? 
2E 



L P ^ \ r 2 - n\ 



(13) 



rf(lnp) 
dr 



-1 



where r\ is the radial position of the L (77) resonance; and r 2 is that the L (77) one, and L p 
6 a p is the mixing angle between the two flavors a and (3 (a,[3=e,fi,r). This condition is satisfied for 
both L-L and 77-77. 



Neutrino Fluxes from Supernovae 

Having constructed the level crossing for the various scenarios, we can now extract the final z^-flux in 
terms of Fi (or F^) in function of the original flux. We can then rewrite it in terms of the flavored 
i/-flux using the identity: 

Fi = Ei Wu\ 2 Fi and F T = EiWuf F- t (14) 

where l=e,fi,r; i=l,2,3; and Uu are the mixing matrix elements. The final flux is thus given by the 
following general relation: 

Fi = Ei auF? + a iT F-° 

F^EianFf + ^Ff 

6 We will denote them as II- a and Il-b. 
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Fig 3: The level crossing diagram for the case of normal mass hierarchy within LMA scheme. 



Fig 4: The level crossing diagram for the case of normal mass hierarchy within LOW scheme. 
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Fig 5: The level crossing diagram for the case of inverted mass hierarchy within LMA scheme. 



Fig 6: The level crossing diagram for the case of inverted mass hierarchy within LOW scheme. 
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where an is the probability that a fi, which is produced inside SN core, leaves the supernova as a 
vacuum eigenstate Vi. Inserting Eq. ()15|) in Eq. I|14)). we find for the flavored neutrino flux: 




Pee Pee 1 Pee Pee 

Pee Pee 1 Pee Pee 

\ I- p ee - Pee 1 ~ Pee "Pee 2 + p ee + p e g + p Ee + Pee 



( F° e 

Fi 

\ K 



(16) 



where p ee is the survival probability of v e , i-e, the probability that a v e doesn't change during the core 
collapse; p^e the probability that a v e leaves SN as a z> e ; the index (°) refers to the original fluxes. It 
is those various matrix elements which characterize the probability of SN neutrino conversion. Those 
probabilities satisfy the condition: 



EaPc/3 = J2f3Paf3 = 1 



(17) 



which are functions of the probabilities aus and U e iS, and a, (3 = e, e, x. 
Case of Normal Mass Hierarchy 

For the LMA scheme: We find from the level crossing for the case of the normal mass hierarchy 
(Fig. 4), that the final flux is given by 



F 1 = (1 - p L )PHFi + (1 - (1 - Pz)p B )F2, F 2 = (1 - PS )F§ + p H F° 
F 3 = F°, F l = p RPL F§ + (1 - PhPl)F2, F- 2 = F°, F 3 = F° x 

then the SN ^-conversion probabilities are given by 

Pee^W- 3 , Pee' ^0.65 X (1 - p L ) p H + 0.35 (1 - PR ) 
Pee ^ 0, p S e ^ 0.65 X p L p R 

For the LOW scheme: According to Fig. 5, the final flux with the LOW scheme is given by 

Fi = P R F° + (1 - p H ) Fi, F 2 = F£, F 3 = p L F° + (1 - p L ) F° 
Fl = PffF§ + (1 - Ph) Fg, F- 2 = (1 - p L ) F° e + p L F°, F 3 = F° 



(18) 



(19) 



(20) 



and therefore the SN ^-conversion probabilities are 



Pee ~ 10 3 X PL , pee ~ 0. 64 X (1 - p R ) 
Pee ^ 0.36 X (1 - pi) , pee ^ 0.64 X p R 



(21) 



Case of Inverted Mass Hierarchy 

For the LMA scheme: The final flux is given in this case (See Fig. 6) by: 

Fi = Pn F£ + (1 - Pr) Fi, F 2 = F°, F 3 = p L F° + (1 - p L ) F° 
Fi = P L F° X + (1 - pi) Fi, F 2 = (l-p H ) F° x + p H Fi, F 3 = F° x 

the SN ^-conversion probabilities are given by 

Pee Ot 10~ 3 X p L , pee ^ 0.65 X (1 - p R ) 
Pee ^ 0.65 X (1 - p L ) , p^ ^ 0.35 X P g 



(22) 



(23) 
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For the LOW scheme: We find for the final flux from Fig. 7: 



F, = p L F° + (1 - p L ) F2, F 2 = I*, F 3 = (1 - pg) F9 + pgF° x 
F 1 =p z F£ + (l-p L )FS, F- 2 = (l- PR )F- + PR Fi, Fg = F% 

and the SN ^-conversion probabilities are given by: 

p ee ~ 0.64 X p L , ^ 10~ 3 X (1 - pg) 

Pee ^ 0.64 X (1 - PL ) , p S e ^ 0.36 X pg 

Note that any significant value for p e g or p Se , is a strong signature of the spin flavor precession effect 
on the SN neutrino dynamics. 



4 Neutrino Conversion Effects on the Mass Spectrum 

Let us discuss the signatures of the various scenarios as far as the emerging neutrino signal, as well as 
their discriminative power at the future neutrino detectors. 

For future supernova neutrino burst, the present detectors can gives, for a typical supernova, at Super- 
Kamiokande 2 about 5000 v e events, and a few hundred events can be detected in both SNO 3 , LVD 
|25j and MACRO 4 . In this paper, we are interested on the features of the final neutrino spectra that 
are required for the identification of the neutrino mass spectra. The effects of neutrino mixing and 
magnetic moment on the final neutrino spectra can be observed through: a) the partial or complete 
disappearance/appearance of the v e neutronization peak, b) the appearance of soft, hard or composite 
spectra of v e and v e ] and c) the Earth matter effects on both u e and v e spectra. Let us estimate these 
effects on the observed spectra at the Earth detectors. 

a) Neutronization Peak: It comes at the first stage of core collapse and corresponds to a and the 
observed signal during the of the neutrino burst first few milliseconds duration. In the absence of 
neutrino conversion, the dominant signal are z^ e s, are produced by the electron capture on protons and 
nuclei while the shock wave passes through the neutrinosphere |14| . Since the original flux is made 
of v e -, the final observed fluxes give a direct measurement of the extent of conversion of v e into the 
other neutrino species. It is thus clear that the neutronization peak is proportional to the i^ e s survival 
probability. 

b) The Nature of the Final Spectra: due to the difference between the interactions strengths of 
the various neutrino species with matter, their average energies differ and are given by |26| : 

(E° e ) = 10 ~ 12MeV, (E? e ) = 14 ~ 17MeV, (E° x ) = 24 ~ 27 MeV (26) 

E means the original average energy (i.e. in the absence ^-conversion). Thus finding for example 
{E Ue ) > (E„ e ), would be is a signature of neutrino conversion, since it implies that the contribution of 
the converted original hard v x spectrum to the final v e flux is significantly larger than its contribution 
to the final v e flux. On the other hand there is another effect which has to be accounted for: the 
fact that the neutrino interaction cross section increases with energy, the neutrino spectra from the 



14 



cooling stage won't be exactly thermal, but will get pinched. We can account easily for this pinching 
effect by parametrizing the original spectra as HZ] : 

^> %W?/r.-,,) (27> 

where T{ and % are given by 

T e ps 3 ~ 4MeV , T F ss 5 ~ 6MeV , T x 7 ~ 9MeV 

The final z^ e (f/ e ) spectrum can be qualitatively divided into three types: 

1) The original "soft" spectrum of the v e (v e ) (corresponding to the survival probability p ee = 1 

(Pee = l))- 

2) The "hard" spectrum of v x (corresponding to the survival probability p ee = p eg = (p S e = Pee = 0) 
which would be the case when there is a complete interchange of spectra, i.e. v e <-> v x {v e *-* u x ); and 
not v e <-> v e ). 

3) The "composite" spectrum, which is a mixture of the original soft and the hard spectra is comparable 
in proportions. There are other cases which are difficult to distinguish like the case where the v e (v e ) 
spectrum contains only v e and u x , v e and u x , or all species. In order to distinguish among them, we 
denote the spectrum containing: v e and v x by compo-1, the one containing v e and v x by compo-2; 
and the one with all the species by compo-3. In some cases, the appearance of one of these case is 
a strong signature of the spin flavor precession effect, namely when the ^ e -spectrum contains v e and 
vice- versa. 

c) the Earth Matter Effect: 

The SN neutrinos, in order to reach the detector, have to go through various amount of the Earth 
material. This amount depends on the direction of the SN with respect to the Earth as well as the 
time of the day. This effect which could modify significantly the flavor composition of the flux may 
be sought by comparing the signal of a future SN from detectors at different geographical locations. 
Certain features of the energy spectra could also reveal this effect even from the observations at one 
detector. 

Neutrinos are expected to be arriving at the surface of the Earth as mass eigenstates (vacuum eigen- 
states); where they oscillate in the Earth matter, i.e, they lose their coherence. The possibility of 
the Earth matter effect's observation depends strongly on the differences between the an and also the 
differences between the Pj e parameters (for more details, see Appendix C). 
Let us now summarize the results within all previous scenarios in the following table: 



5 Conclusion 

We attempted in this paper to check the sensibility of the supernovae final flux of neutrinos to both 
the magnetic field and matter effect assuming three active neutrinos. We have taken the value of \U e 3\ 
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Region 

Of llB 1 n 


SNP solution 


Neutro- 
Peak 


Spectrum of 

l> e v e 


Earth effects 




normal 

LMA 


v x 


hard 


compo-2 


X 


/ 


I 


inverted 


v x 


hard 


compo-2 


x 


/ 


normal 

LOW 

inverted 


v x 


hard 


compo-2 


X 


/ 




" e 5 " x 


compo-1 


compo-2 


/ 


/ 




normal 

-a LMA 

inverted 


v e , v x 


compo-2 
compo-2 


hard 
compo-2 


/ 
/ 


/ 
/ 


II 


normal 

-b LMA 

inverted 


Ve-,V X 

v e , v x 


compo-2 
com'no-2 


compo-2 
corrvno-3 


/ 
/ 


/ 
/ 




normal 

LOW 

inverted 


V e , v x 


compo-2 


compo-2 


/ 


/ 




v e ; " x 


compo-1 


compo-2 


/ 


/ 




normal 

LMA 


v e , v x 


compo-2 


hard 


/ 


X 


III 


inverted 


v e , u x 


compo-2 


compo-2 


/ 


/ 


normal 

LOW 

inverted 


v e , v x 


compo-2 


hard 


/ 


/ 




v e , v x 


compo-1 


compo-1 


/ 


X 



Table 1: The dependance of the final spectra on: 1) the SNP solution, 2) the mass hierarchy and 3) the range 
of /jtB± Q . "soft" in the column refers to the original v e (u e ) spectrum and "hard" refers to the original v x 
spectrum. In the Earth matter effect column, / and x indicate the possibility of significant effects. 

to be just below the experimental upper bound given by CHOOZ [2J|) which makes the transition at 
the H layer completely adiabatic. We then divided the whole range of values considered for fJ,B± 
into three regions /, //, and I/I according to the adiabaticity at the H layer. We could then make 
some predictions on the conversion effects for supernova neutrinos. The predictions differ for the 
different schemes, which opens up the possibility of discriminating between them using data from 
future neutrino bursts. We studied the possibility of observing the conversion effects through: (i) The 
partial or complete change of the flavor of the neutronization peak, (ii) The appearance of a hard or 
composite v e and/or v e spectra due to the conversion effect instead of their original soft spectra. (Hi) 
the Earth matter effect from the conversion in the Earth material. We found that indeed neutrino 
conversion does change significantly the spectrum shape. Let us now summarize the salient features 
of the final flux as it appears from the last table: 

1. The appearance of a hard v e makes the mass hierarchy to be normal. 

2. A hard v e spectrum makes /J,B ± to be in the first region, while a hard v e spectrum exclude the 
first region. 

3. The absence of the Earth matter effect in the v e channel implies that \xB ^ to be in the first region; 
while the absence in v e one implies that it must be in the third region. 

4. If v e contains all spieces (i.e. is e , v e and v x ), the SNP solution must be LMA, the mass hierarchy 
is inverted and [iB± lies between 1700fj> B G and 2400fj> B G. 

5. If both spectra of v e and v e contain only v e and u x , the mass hierarchy will be inverted, the SNP 
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solution will be LOW, and (xB±_ will be larger than 2400/i B G. 

Notice that some possible observations can rule out all of the various scenarios above, like the appear- 
ance of a soft v e (or v e ) signal, which would exclude all the above scenarios, since then the value of 
\U e z\ 2 would be less than 10~ 3 . 

The final neutrino spectra can thus help in resolving three main kinds of ambiguities that remain 
to be resolved with the current data: (i) the solution of the SNP, (ii) the type of mass hierarchy 
(sign of Am^), and (to) probe the magnitude of the [iB j_ value, assuming \U e3 \ 2 to be 10~ 3 . The 
implications of the results of this work will depend on when will the next neutrino burst from a 
Galactic supernova be detected. On the other hand, there is a good chance that within the next few 
years the present, ongoing, or future planned experiments will allow us to identify the specific solution 
of the solar neutrino problem considered in this work. This will significantly diminish the number 
of possible schemes and will allow us to further sharpen the predictions of the effects for supernova 
neutrinos. 
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A In which Channel does the L Resonance Occur? 

Taking the Hamiltonian in the (v e ,is a ) basis: 

( W K \ 

(29) 



1 W K 
2E\ K Qz 



where v a is any combination of the non-electronic neutrinos, K and W are given in function of the 
various Am 2 's; and Q is a function of fl^'s and Y e . We know that the resonance occurs at z res = 
if z res < 0, then the transition occurs between {p e , u a ) instead of {y e ,u a ). In general Q is always 
positive, thus we should check the sign of W. 

W represents the eigenvalues of the Hamiltonian submatrix, which corresponds to the non-electronic 
states {y^^v-j). (see. Sec II); In order to probe the sign of these two values in the general case, we take 
a 2 x 2 symmetric matrix, similar to the mass matrix in the subspace (z/ M ,i/ T ), as follow 



e A 
A a 



(30) 



where e, A and a are reels; it's eigenvalues are: \{e + a ± \j(e — a) 2 + 4A 2 ); then one can see that the 
two eigenvalues have: 

• the same sign if |e<j| > A 2 . 

• different signs if \ea\ < A 2 . 

The small value, in absolute value, leads the resonance at the L layer, while the largest one leads to 
the resonance at the H layer. Thus, we are led to look at the sign of the smallest value, replacing e, 
A and a by their values in our case, and neglecting terms multiplied by ^^i 1 , then one can find the 
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condition that the resonances at both L and H layers, occurs in the same channel, it can be written 

as 

sin 2 fl 13 <^cos2# (31) 

this condition is satisfied in the case of inverted mass hierarchy with both LMA and LOW scenarios. 
But for the case of normal mass hierarchy, it is satisfied only within LMA scheme. 



B The Adiabaticity of Neutrino Conversion 

For neutrinos traveling through a non-constant density medium, it is proved that the adiabaticity is 
always satisfied except for the cases around the resonance layer, where it must be studied more thor- 
oughly. The resonance, in general, can be characterized by the matter density value p, and therefore 
the distance r from the SN center, when the difference between the two Hamiltonian eigenvalues, 



Mf M\ 
IE ~ ~2E 



, is minimal, i.e. 

d 



w(A#e//| res = (32) 

where / is either the density p or the travelled distance r. Considering the 2-v scheme (Like v e ,v^) 
case for the MSW effect, H can be written as: 

IE \ m\ ) 2 1 \ Y e -1 ) ' 
and the mixing matrix U, is given by 

U=( COS6 S[n6 ) (34) 
y — sin 6 cos 9 J 

Neglecting the terms proportional to the unity matrix, we find 

/00\ + /l0\ / ap + Asin 2 (9 #sin26>\ 
H = U [o A)' J, + aP { oj = ( |sin2, Los'e) (35> 

where A = 2 2E 1 . Then the difference between the two eigenvalues is: 

\J {cap - A cos 26f + A 2 sin 2 26 (36) 

which takes its minimum at 

P=#cos2# (37) 
This corresponds to the equality of the diagonal elements. 

If we try to redo the same operation but writing this time H in the vacuum state basis (^1,^2), we 
find 

rr rrirrrr ( \ ( COS 2 6 \ Sm 26 \ . , 

H = U ] HU = + ap \ i 2 . (38) 

I A J H \ isin20 sin 2 ^ J V ' 
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The difference between the eigenvalues takes its minimum at the same value of p, given in Eq. Q37JI . 
because they are different representations of the same operator H, which doesn't correspond to the 
diagonal elements equality 

P = ^kre (39) 
We will obtain the same result if we write H in any other basis, i.e., the minimum of AH e ff doesn't 
correspond the diagonal elements equality. 

One can deduce that the minimum of AH e ff corresponds to the diagonal elements equality only when 
the off-diagonal elements are constant, i.e. /^-independent which is realized only in the flavored basis. 
There are some cases, like the one in the presence of magnetic moment interaction, where the flavored 
basis is not the best basis to deduce the resonance density directly, i.e. by equating the diagonal 
elements. 

Let us consider our case, from Eq. the minimum of AH e ff is given by 



£,y/(b-CZy + 48*Z* (40) 

which corresponds to 



fee 



^res — c 2 +4s 2 

but we take it -; since s <C c, then one can write 



(41) 



c 



l-7.46xl0" 14 x (i^f) ) (42) 



2 



thus the layer radius r res is corrected by Ar res , where Aires _ _4_21 x 10 5 x ( J 3 , then the 
adiabaticity parameter will be corrected by the factor 

3 v ( 



£ = 1 - 1.24 X 10~ 13 X (43) 



and the jumping probability becomes 



^(corrected) = ^ ^ 



suppose that Pf =0.4, in order to be corrected by 0.1, i.e. p( C0Trected ) =0.5, we need that the quantity 
fiB± to be larger than 1.4 x 10® fi B G; which is a huge value non-available in the supernovae. We 
conclude that this correction doesn't affect the neutrino conversion in our work, and we have thus 
neglected it. 

C The Earth Matter Effect 7 

If neutrinos reach the Earth detectors without interacting with matter, the signal at the detectors is 

given by Eq. (|16[). Since they reach detectors through the Earth matter, the flux given by Eq. (fT6|) 

will be modified in general. 

7 We are here closely following the treatment of adopting it to our case by taking into account the possibility of 
the v — v transition. This effect is studied with more details for difference scenarios in |28|. 
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Let Pi e (Pie) be the probability that a vacuum mass eigenstate V{ (i/j) entering the Earth reaches the 
detector as a u e ipej- The flux of v e (v e ) at the detector is 

Fe D = Ei^eFi, Fe D = E<P^ (45) 

Inserting Fi, we get 

Pe — F e X^i Pie&ie ~\~ Pg X;i Pie&ie ~\~ Px X]j Pie (1 &ie &i 



ne ^le ) 

Pe — Pe X^i PieO'ie ~i~ Pg X]j Pie®ie Px X;i Pie (1 ^Tey 



(46) 



where the a's are defined in Sec. 3. Then one can write 

Pee = X;i Pie a iei Pee = X^i Pie a iei Pee = X^i Pie a iei Pee = X^i Pie a ie (47) 



Similarly we have 



which enables us to write: 



F e = VeeF° e + PeeF§ + (1 - Pee ~ Pee) F° 
F-e = p-eeF? + p- ee F° + (1 - Pee ~ Pee) F° x 



rpD pD po , pDpp , / 1 _ pD _ p£A po 

x e x ee e 1 x ee e 1 v ee ^ ee) x 

pD _ pD po I pD _ p£> _ pD\ po 

x e - t ee i e 1 - 1 ee e 1 \ ee x ee/ 5 



(48) 



(49) 



From those expressions we deduce: 



F e D - F° = {F° - F£) (P e D e - Pee) + (F§ ~ F° x ) (Pg - Pee) 

F e D -F£ = {Fi - F£) {Pg - p-e-e) + (F e ° - F°) [Pg - Pse ) 



(50) 



so the Earth matter effect can be quantified by the various differences (Pg—p ee ), (Pg—pee), (Pg—Pee) 
and (Pg — Pee), which equal; 

Pee Pee — X] &ie(Pie |^ei| )> Pee Pee — X] ®ie(Pie |^ei| ) (51) 
Pee Pee — X^ &ie(Pie |^ei| )i Pee Pee — X^ &ie(Pie |^ei| ) 

One can finally write Pg — p ee as 

Pg - pee = a le (Ple ~ \U e l\ 2 ) + a 2e (P 2 e ~ \Ue2?) + a 3e (P 3e " \U e ^) (52) 

We obtain similar expressions for the other difference terms. We notice that the last term in Eq. (|52|) 
is negligeable due to the very small depth oscillation of inside the Earth 

P 3 e-|^e3| 2 < 10- 3 (53) 

Taking into account the relations X] |£4i| 2 = 1 and X^j Pie = Yl Pie = 1; one ca n write 



Pg - Pee = (a-2e ~ CUe) (P 2e ~ | ^e2 1 2 ) , Pg ~ Pee = («2e ~ 0>le) (^2e ~ We2\ 2 ) 
Pg ~ Pee = (a 5 g - a- ls ) (P^e ~ \ U e2\ 2 ), Pg ~ Pee = («2 e ~ Oje) ( P 2e ~ We2\ 2 ) 



(54) 
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finally, the differences Fj? - F° and Ff - F? can be written as 

F e D - F° « [{F° e - F°) (a 2e - a le ) + (F? - F°) (a 2e - - a le -)] (P 2e - \U e2 \ 2 
Ff-Fi* [{F9 - F°) (a 2g - a Ig ) + (F e ° - F°) (a 2e - a Ie )] (P 2g - \U e2 



2 (55) 



In general, when the signal from two detectors D\ and D 2 are compared, we get the following flux 
differences: 

F Dl _ F D 2 „ {i po _ po) (a2e _ aie) + {F o _ po^ (a2g _ aig)] ( p(l) _ pg)) 

if 1 - if 2 ^ - F, ) (a 2e - - a Ig ) + (F e ° - F°) (a 2e - o Ie )] - p£>) 

According to this relation, the difference F® 1 — F® 2 (F B 1 — F® 2 ) is factorized: it is proportional 
to the difference of the Earth oscillation probability P 2e (P 2S ) at the two detectors and a function of 
supernovae oscillation probabilities a's and the difference in the original fluxes of v e (v e ) and v x . 
Let us consider these factors separately: 

1 ) P 2 * — P^f and P^_ — P^- : if the neutrino trajectory go through only the Earth's mantle, one 
can use a constant density approximation, which gives 

He ~ ?2e = sm20r 2 sin2(0f 2 - 0i 2 ) 



Pff - pg } = sin20gsin2(0g - 9 e2 ) 



sin 2 ( *dA _ gin 2 f 
sin 2 f ^1 ) _ sin 2 ( W 



(57) 



2e - 1 2e ~~ ^^ a ^^ a 

here # m and L m are the mixing angle and the oscillation length inside the Earth respectively, and di 
is the distance travelled by neutrinos inside the Earth before reaching the detector Di. 

2) The second factors is a summing of two-factor terms. Let us try to see how they behave: 

* (F° — F°) and (F g — F°) : Since the v e (v e ) spectrum is softer than the u x spectrum, and the 
luminosities of both the spectra are similar in magnitude [2JJ], the term (F° — F°), (and therefore 
(Fg — F°)), is positive at LOW energies and becomes negative at higher energies where the u x flux 
overwhelms the v e {p e ) flux. Therefore, the Earth effect has a different sign for LOW and high 
energies, and there exists a critical energy E c (E c ), such that F°(E C )= F°(E C ), (F?(E C )= F°(E C )), 
where this change of sign takes place. Since the cross section of the neutrino interactions increases 
with energy, the Earth effect is expected to be more significant at higher energies (if all the other 
factors are only weakly sensitive to the neutrino energy). 

* (fl2e — «ie) ; («2e ~~ a ie)> ( a 2e — ais) or (a 2g — aj g ): They represent the differences between the an 
probabilities and can be extracted from Eq. (|18j). Eq. ((20)) , Eq. ((221) and Eq. (|24[). They are seen to 
be dependent on the mass hierarchy as well as the specific solar neutrino solution (LMA or LOW). 
Non vanishing values for this difference make the observation of the Earth matter effect possible. 
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